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Limits of complete equilibration of fragments produced in central 
Au on Au collisions at intermediate energies 
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Abstract. Experimental data related to fragment production in central Au on Au collisions were analysed 
in the framework of a modified statistical model which considers cluster production both prior and at the 
equilibrated stage. The analysis provides limits to the number of nucleons and to the temperature of the 
equilibrated source. The rather moderate temperatures obtained from experimental double-yield ratios of 
d,t,^}ie and ^He are in agreement with the model calculations. A phenomenological relation was established 
between the collective flow and the chemical temperature in these reactions. It is shown that dynamical 
mechanisms of fragment production, e.g. coalescence, dominate at high energies. It is demonstrated that 
coalescence may be consistent with chemical equilibrium between the produced fragments. The different 
meaning of chemical and kinetic temperatures is discussed. 

PACS. 25.70.Mn , 25.70 Pq ,25.75.-q 



1 Introduction 



The mechanism of fragment production at intermediate 
energy nucleus-nucleus collisions is a long standing prob- 
lem which is important for many investigations, such as 
the study of the nuclear liquid-gas phase transition, the 
behaviour of nuclear matter under extreme conditions, 
and, in general, the nuclear equation of state. In this paper 
we concentrate on central collisions which deposite a large 
amount of energy into nuclei, and produce a fast explosion 
of nuclei into many fragments. Presently, there is evidence 
that at low projectile energies of '^10-100 A-MeV the 
fragment production via multifragmentation of thermal- 
like sources is the dominating process 1,2,3,4 . With in- 
creasing energy dynamical effects, such as collective flow, 
becomes prominent. This implies that the fragmentation 
mechanism changes from a statistical to a dynamical one. 
However, the description of intermediate mass fragment 
(IMF) production as a result of the break-up of an equi- 
librated source with collective flow is very successful, as 
shown in many publications We believe that such 

good descriptions were possible since the hypothesis of 
chemical equilibrium between different kinds of fragments 
is adequate for these reactions. 

In this paper we analyse data from 100 to 1000 A-MeV. 
Some of these data were already analysed with different 
dynamical and statistical approaches. The present analy- 
sis involves new degrees of freedom. In particular, it em- 
phasizes the importance of isospin characteristics of pro- 



duced fragments for the determination of the reaction 
mechanism. 



As suggested by many theoretical and experimental 
studies P||5|[7| | 5| PllUlllH the fragments may be produced 
in a fast initial ('preequilibrium') process as well as at 
the full equilibration stage. Since the number of nucle- 
ons participating in the thermal-like source can decrease 
rapidly with the beam energy, dynamical processes of frag- 
ment formation should contribute essentially at higher 
energies. The mechanism of coalescence of nucleons into 
fragments is applied here for a complementary descrip- 
tion of light fragment production to the statistical ap- 
proach. In the following we discuss a relation between 
coalescence and statistical approaches and we point out 
that the coalescence mechanism may simulate the chemi- 
cal equilibrium conditions. In this respect, we pay special 
attention to the chemical temperature and study its cor- 
relation with the light charged particles (LCP) velocities. 
However, our main aim is to determine properties of a 
completely equilibrated source, which is mainly responsi- 
ble for IMF production. We emphasize, that the conclu- 
sions about equilibrium and nonequilibrium contributions 
concerns the production of fragments from nucleons only. 
We do not consider the problem of the nucleon thermal- 
ization and refer to experimentally selected central events 
which belong to an ensemble of nucleons in some degree 
of equilibration. The knowledge of the relation between 
equilibrium and non-equilibrium mechanisms of fragment 
production in dynamical processes is important for many 
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applications both in nuclear physics and astrophysics. For 
example, similar processes are expected during the fast 
synthesis of elements in the early universe and in super- 
nova explosions |12j . 

2 Chemical temperatures evaluated from 
fragment data obtained in central Au+Au 
collisions 

Data obtained by the FOPI collaboration at 100, 150, 
250 and 400 A-MeV pitl^ and data from the EOS col- 
laboration |16II17II18| at 1 A-GeV are analysed. Central 
events were selected by the criterion ERAT as described 
in ref. |E| whereas the event centrality of the 1 A-GeV 
data was determined by multiplicity cuts (see ref. |17j). 

It is commonly accepted that the kinetic energies of 
fragments can be represented as a sum of thermal and 
collective components. The collective motion (radial flow) 
is an important ingredient of Au on Au collisions at in- 
termediate energies and influences strongly the fragment 
energies. Kinetic energy distributions for central events 
have been analysed in the framework of the 'blast model' 
which is described in detail in ref. 13j. Here, we recall 
briefly only the main aspects. The collective energy stored 
into radial flow was determined by velocity profiles and an 
ansatz for the velocity distribution ^] of the fragments. 
The corresponding kinetic energy distributions were re- 
produced if the collective energy ScoU amounts to 62±8 % 
of the center-of-mass energy Ec.m. available in the colli- 
sion. Then, the energy conservation requires that 

£coU + £th — Ec.M. + Q (1) 

where Q is the Q-value of the reaction and £th is the ther- 
mal energy. In the limit of classical statistics the temper- 
ature is determined by the multiplicities N of the emitted 
particles using the non-relativistic expression 

etu^\-{N-l)-T*. (2) 

The temperatures found by this approach are T*= 
17.2±3.4, 26.2±5.1 and 36.7±7.5 MeV for the beam en- 
ergies 150, 250 and 400 A-MeV, respectively. Below we 
label them as 'kinetic' temperatures Tkin- It was empha- 
sized that the temperatures T* are 'effective' in the sense 
that they are not the temperature at freeze-out time, since 
the observed multiplicity may be raised due to late par- 
ticle decays. In ref.^Hl, these temperatures were used as 
input for the statistical multifragmentation models QSM 
pn) . WIX in] and SMM fF3^ by assuming that the whole 
entire mass of the colliding Au nuclei undergoes thermal- 
ization. This assumption fails to reproduce the fragment 
multiplicities. The abundances of heavy clusters are un- 
derestimated up to 4 orders of magnitude (e.g. for Z=8 
fragments at 400 A-MeV.) In ref.^^l it was discussed that 
variations of the radial flow energy, the freeze-out densi- 
ties and level density parameters within reasonable physi- 
cal limits cannot account for such large deviations. Hence, 



the question arises about the applicability of established 
statistical models to cluster production in the midrapidity 
source and the meaning of the nuclear temperatures men- 
tioned above. The present paper is intended to help to 
disentangle this 'puzzle'. First we ask whether these tem- 
peratures can be confirmed by other data available for the 
same collision system Au+Au measured with the same ap- 
paratus. Such possibility offer the yields of hydrogen and 
hehum isotopes ^S] which can be treated by the isotope 
thermometry. 

In ref. 23.^ it was shown that in an equilibrated sys- 
tem the double-yield ratio {R1/R2) of isotopes is directly 
related to the temperature of the corresponding grand 
canonical ensemble: 

^ ln(a.(i?i/i?2))- 

For the consideration of hydrogen {d, t) and helium iso- 
topes {^Hef^ He) one needs to fix the parameters 6=14.32 
MeV and a=1.59 ||24j which include the binding ener- 
gies, masses and spin degeneracy. This isotope thermome- 
ter has proved to be successful in many applications. In 
particular, for the first time it was possible to establish 
experimentally the nuclear caloric curve |25 . 

Here, we refer to the isotopic yield ratios and kinetic 
energy distributions of d, i, ^He and '^He measured with 
AE/E telescopes within the CM. polar angle range of 
60° < Oc.M. < 90° Most of the projectile fragments 
are expected to be suppressed within this angular cov- 
erage. The corresponding data sets at 100, 150 and 250 
A-MeV were obtained from central event samples selected 
by the criterion ERAT5 13 . The integration of the C.M.- 
kinetic energy spectra delivered the intensity of deuterons, 
tritons, '^He and ^He from which the ratios Ri and R2 and 
the corresponding isotope temperature Tiso were deter- 
mined. The obtained values (T^so) = 6.34 ±0.50, 7.8 ±0.8 
and 11.51±1.58A/el/for 100, 150 and 250 A-MeV, respec- 
tively, are about three times smaller than the correspond- 
ing kinetic temperatures. More details of the distribution 
of Tiso in the freeze-out volume are required. Here, we de- 
termined from the LCP spectra presented in ref. the 
isotope temperature in dependence on the velocity. 

The CM. kinetic energies Ekin were transformed into 
the particle CM. velocities using the relativistic relation 

v_ _ \/ EkinjEkin + 2mc^) 
c Ekin + mc^ 

where m is the corresponding LCP mass, c is the light 
velocity. The spectra given in ref. have an equidis- 
tant energy binning for all particles. But, the velocity di- 
visions become varying for different particle masses after 
the abscissa transformation. In order to cover the same 
range of velocities we analysed the kinetic energy distri- 
butions to Efc,„(d) < 110 MeV, Ektni.t,^ He) < 150 MeV 
and Efci„(^iJe) < 200 MeV for incident energies 100 and 
150 A-MeV. The corresponding limits are 170, 250 and 300 
MeV at 250 A-MeV. Equal reference velocities v,; for the 
four LCP's were found by an appropriate interpolation 
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of the velocity distributions derived within these limits. 
The maximum velocity up to which yields of all LCP's 
were available was around 0.3-c at 250 A-MeV and 0.6-c 
at lA-GeV, respectively. Then, the isotope temperature 
Tiso was recalculated for subsequent velocities v^ by means 
of equation (3) with the ratios 

Ri^Y{d,Vi)/Yit,v^) 
R2 = Y{^He,v,)/YCHe,v,). (5) 

The mean temperatures evaluated within the limited 
velocity ranges are (Ti,o) = 5.30, 6.10 and 9.14±0.60 MeV, 
respectively, for 100, 150 and 250 A-MeV. They are smaller 
than the values derived from the yields integrated over the 
complete energy spectrum, since the contributions of frag- 
ments with the highest energies are missing. However, this 
treatment is self-consistent for this subset of events, where 
the collective flow dominates over Coulomb and thermal 
energies and determines mostly the fragment velocities. 

As shown in Fig. 1 the temperatures Tiso become larger 
with increasing particle velocity, however, they are below 
the kinetic temperatures T* of ref. ^HI- This finding sug- 
gests that the obtained data cannot be described by one 
temperature characterizing the complete equilibrium. An 
explanation of this phenomenon could be that fast (pree- 
quilibrium) nucleons are emitted earlier at very high tem- 
peratures and carry away excess energy, while slow nu- 
cleons form an equilibrium- like source The observed 
light fragments are produced during the cooling process 
and the distribution of the chemical temperature versus 
the fragment velocity reflects evolution toward equilib- 
rium. Figure 1 shows that in all cases the temperature 
Tiso has some saturation at low velocities. This could be 
an evidence for reaching thermalization of nuclear matter 
at these velocities. Obviously, this thermal source can be 
characterized by very moderate chemical temperatures. 

It is a general observation (e.g., see that in in- 
termediate energy collisions the kinetic energy spectra of 
fragments cannot be reproduced by using an unique tem- 
perature. In the following we propose a model to describe 
these feature. 



3 Description of the model 

An appropriate way to describe processes involving many 
particles is the statistical approach. The system charac- 
terized in the initial stage by nonequilibrium distribution 
functions evolves towards equilibration as a result of many 
interactions between the particles. In this process the sys- 
tem runs through different states. The first one can be 
considered as equilibration of the one-particle degrees of 
freedom. The following evolution toward total thermaliza- 
tion can be considered as involving of higher order parti- 
cle correlations. For finite expanding systems the degree 
of equilibration depends on the reaction type. It is ex- 
pected that the equilibration is less effective for nucleons 
at the surface than for nucleons deep inside the freeze- 
out volume. As a result, the mechanisms of the fragment 
production may be different also. 



lOOAMeV 150AMeV 250 AMeV 




5 100 5 100 5 10 



velocity (cmns"^) 



Fig. 1. Isotopic temperatures evaluated from the data in 
ref. 1151 as a function of the CM. velocities of d,t,^He and 
'^He. Dots: data. Lines: calculations with the code SMMFC 
including contributions from the fast stage (coalescence). 



In the following we are going to interpret the experi- 
mental data within the framework of the expanded version 
of the statistical multifragmentation model 1 , 22 which 
treats also LCP's emitted prior to the equilibration. The 
model phenomenologically includes collective motion (ra- 
dial flow {F)) and, optionally, composite particle produc- 
tion by a new coalescence (C) algorithm. ^ The subsequent 
application of these mechanisms is aimed at simulating 
the most important physical processes in many-particle 
systems. 

The model parameters have to be derived from fits 
to the experimental data as described in section 4. It is 
important to use as many as possible observables for the 
analysis, and to ensure that the deduced parameters and 
their dependence on the beam energy are qualitatively 
consistent with general predictions of dynamical calcula- 
tions 7 . 

Total energy and momentum balance is used for the 
implemented processes (except 7-emission). First we con- 
sider the quantities being constrained. A projectile with 
mass number Ai and charge Zi collides at beam energy 
Ebeam with a target {A2, Z2) resulting in the center-of- 
mass energy Ec.m. available for the total system with 
Ao=Ai+A2 and Zq~Zi+Z2. Since the equilibrated stage 
is rather well understood, we start the analysis of exper- 
imental data related to IMF's which are assumed to be 
produced only in the equilibrated (i.e. thermal) source. 
This source is parametrized by (i) the mass number As — 
Arei ■ Aq (where Arei is the relative source size) and corre- 
sponding charge Zs = As-{Zo/Ao), (u) the thermal excita- 
tion energy E* and (iii) the collective energy per nucleon 
Efiow In the analysis of the 1 A-GeV data the energy re- 
leased due to pion production was taken into account. Pjr 
is the part of Ec.m. taken away by pions. The remaining 
matter {Apre — Aq — As, Zpre = Zq — Zs) is assumed 
to be carried away by fast nucleons and LCP preequi- 
librium emission. In this case the change of the binding 
energy is given by the corresponding values for projec- 



^ henceforth the sign SMMFC denotes this code version 
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tile, target, thermal and preequilibrium sources: AB 
Bi + B2 — Bg — B 



pre ■ 



From the conservation of the to- 
tal energy follows the energy available for preequilibrium 

Ptt ■ Er:_ M. — E* — Efiow + AB. 



emission 



E„ 



E, 



CM. 



3.1 LCP emission at the fast stage. 

Since the cooling during the expansion process is very fast, 
nucleons have no time to feel the part of the phase space 
corresponding to the composite particle production. A dis- 
tribution of nucleons in the phase space at some 'freeze- 
out' time is considered as start configuration. Generally, 
any distribution of nucleons in momentum and coordinate 
space after an initial dynamical process is conceivable. 
But, in some experiments, e.g. central nucleus-nucleus col- 
lisions 221 ) it is possible to select samples which are nearly 
isotropic in space and look like thermal events. Therefore, 
in such cases, we can simply assume that the nucleons 
populate the available many-body phase space uniformly, 
i.e. there is equilibration in one-particle degrees of freedom 
without manifestation of collective phenomena. That gives 
rise to a thermal distribution for individual nucleons in the 
thermodynamical limit. In the calculations we consider 
the system characterized by Apre , Zpre and Epre and dis- 
integrate the system into nucleons by taking away about 
7 A-MeV (binding energy). The remaining energy turns 
into the kinetic energies of nucleons which populate the 
whole available many-body momentum space uniformly. 
The procedure developed in [221 is used to generate the 
nucleon momenta. 

A composite particle can be formed from two or more 
nucleons if they are close to each other in the phase space. 
This simple prescription is known as coalescence model. 
Here we use the coalescence in momentum space which 
was recently described and applied in ref. 29 . The basic 
assumption is that a dynamical process, which leads to 
a momentum redistribution, is very fast (nearly instanta- 
neous), so that the coordinates of nucleons are just de- 
fined by their momenta. This is also justified taking into 
account quantum properties of the system since the wave 
functions of nucleons can be described in momentum space 
only. This type of coalescence model has proven success- 
fully by reproducing experimental data (see e.g. [3DII31L 

In the standard formulation of the model it is assumed 
that the fragment density in momentum space is propor- 
tional to the momentum density of nucleons times the 
probability of finding nucleons within a small sphere of 
the coalescence radius pq. For example, in a nonrelativistic 
approximation, from this hypothesis an analytical expres- 
sion can be derived for momentum spectra of coalescent 
clusters: 



numbers A and 1, respectively. This equation disregards 
correlations between different clusters since the conser- 
vation of the nucleon number is not taken into account. 
Therefore, the above formulae is valid only for (A^i) 3> 

{N2) » (iVg) ... 

We developed another formulation of the coalescence 
model. Nucleons can produce a cluster with mass number 
A if their momenta relative to the center-of-mass moment 
of the cluster is less than pq. Accordingly we take |pi — 

Pc.M.I < Po for all i = 1, A, where pc.M. = 3 J2t=i Pi- 
In the following examples the value po ~ 94MeV/c has 
been adopted corresponding to relative velocities Vrei—O-^c 
in agreement with previous analyses |29j . 

In this context we would like to draw attention to a 
problem which is sometimes disregarded in these simula- 
tions. Some nucleons may have such momenta that they 
can belong to different coalescent clusters according to the 
coalescence criterion. In these cases the final decision de- 
pends on the sequence of nucleons within the algorithm. 
To avoid this uncertainty we developed an iterative coales- 
cence procedure. M steps are calculated in the coalescence 
routine with the radius poj which is increased at each step 
j: poj — {j/M.) ■ Po U — l,...,M). Clusters produced at 
earlier steps participate as a whole in the following steps. 
In this case the final clusters not only meet the coalescence 
criterion but also the nucleons have the minimum distance 
in the momentum space. This procedure is unique in the 
limit A4 —>■ 00 and we found that in practical calculations 
it is sufhcient to use A4=5. 

The importance of this coalescence mechanism is demon- 
strated by Fig. 2 which shows a comparison of measured 
isotopic yield ratios d/t and ^He/'^He in the beam en- 
ergy range from 35 A-MeV to 1 A-GeV |15lll6ll27ll28j . 
Especially at higher incident energies the SMM calcula- 
tions without consideration of coalescence underestimate 
strongly the existing data. The so-called ^^He He puz- 
zle' can also be solved by taking into account the coales- 
cence 1211 ■ 

In the calculations we assume that only coalescent par- 
ticles with A < 4 are produced. In principle, one can 
extend the model by considering IMF also, though the 
probabilities of coalescent IMF is small j32ll29j . In this 
case the expected portion of the thermal source will be 
even smaller than with restriction to A < 4. However, this 
paper is aimed at finding an upper limit for the contribu- 
tion of thermal IMF's. As justification of our assumption 
we present a good agreement with experimental data in 
section 5. There are also other experimental features sup- 
porting this assumption. In particular, the maximum of 
IMF production was found at small fragment velocities. 
These velocities correspond to the nearly constant tem- 
peratures Tiso (see Fig. 1) associated with a thermal-like 
source. 



dpi 



An 



A-i 



(6) 



where p„ are the momenta per nucleon. (Na) and (A^i) 
are the mean multiplicities of fragments with the mass 



3.2 Thermal source 

At the freeze-out time of several tens of fm/c there is still 
a lot of relative dense nuclear matter in the center of the 
system in which instabilities (usually associated with the 
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3.3 Relation between coalescence mechanism and 
thermal fragment production. 

An important relationship can be formally established 
between the coalescence and thermal models as far as frag- 
ment production is considered. This will be illustrated by 
a simple statistical case. A system with the total num- 
ber of nucleons Aq disintegrates into fragments with mass 
number A, which are characterized only by their bind- 
ing energy Ba- The fragments are considered as Boltz- 
mann particles moving without interaction in the volume 
V. Then the statistical partition sum can be calculated as 




Pa 



(8) 



Fig. 2. Isotopic yield ratios of LCP in central Au on Au col- 
lisions as a function of the beam energy. Data are displayed 
by symbols: FOPI EOS [T5], MSU [23, Doss |2H1- 

The hatched areas show SMMFC calculations assuming the 
uncertainty of the input parameters given in tablet. Dashed 
lines: Calculations using only the thermal source (without coa- 
lescence) and the same input (without consideration of uncer- 
tainty limits). 



liquid-gas type transition) may occur. At this stage it is 
assumed that thermal equilibrium is reached, which is re- 
sponsible for IMF production. The success of the SMM 
and other statistical models [Tll21ll33| to describe such re- 
actions (see e.g. refs . [T]|2||5| | 3^1 E 1 ) supports this assump- 
tion. 



A microcanonical ensemble of all break-up partitions 
which consist of nucleons and excited fragments of dif- 
ferent masses is considered. If E* and the volume of 
the thermal source are fixed, the statistical weight W of 
a given break-up partition j (i.e. the number of micro- 
scopic states leading to this partition) is determined by 
its entropy Sj: 



Zs)). 



(7) 



The fragments with mass number A > 4 are treated as 
heated nuclear liquid drops but the light fragments with 
A < 4 are considered as massive particles ('nuclear gas') 
having only translational degrees of freedom. The ensem- 
ble of partitions is generated by Monte-Carlo methods ac- 
cording to their statistical weights Wj constrained by the 
conservation laws. The microcanonical temperature Tth 
is found from the energy balance by taking into account 
the Coulomb interaction, binding energies and excitations 
of fragments. After break-up of the system the fragments 
propagate independently in their mutual Coulomb fields 
and undergo secondary decays. The deexcitation of large 
fragments {A > 16) is described by the evaporation-fission 
model, and for smaller fragments by the Fermi break-up 
model EIEII. 



where T is the kinetic temperature characterizing the frag- 
ment translational motion, rn„ « 0.94 GeV is the nu- 
cleon mass and p is the fragment momentum. The prod- 
uct includes all fragments in a partition. The momen- 
tum and the center of mass conservations are disregarded. 
The magnitude Pa is proportional to the probability of 
the formation of fragment A, which can be written as 
Pa = exp{—BA/T) in the case of full equilibration charac- 
terized by the canonical temperature T. However, the pro- 
cess of fragment formation may be complicated, and this 
process may not be related only to the thermal kinetic mo- 
tion. Generally, one can write Pa = exp{—BA/TA), where 
Ta is a parameter ('temperature') related to a given frag- 
ment. One can introduce a Lagrange multiplier fi (like a 
'chemical potential'), which can be found from the con- 
dition '^Ai^A)A = Aq (see ref. J,). Then the partition 
sum can be calculated as 



^-E E n 

Ni=Q Nao=0 a 



{Na)''- 



NaI 



(9) 



Here Na is the multiplicity of fragment A in a partition, 
and 



(Na) = ^A3/2exp 



(10) 



1/2 

with the thermal wavelength At = (27r?i /m„r) . The 
last equation can be rewritten as 

(Na) = {N^ (|^)^"^^/^exp(-|^). (11) 

At this point we can establish a correspondence between 
the thermal and the coalescence models (see also |^). In 
the coalescence model, the fragment multiplicity can be 
determined after integration of equation © . If the nucle- 
ons are assumed to have a Maxwell-Boltzmann distribu- 
tion of the same kinetic temperature one can easily get 



{Na) 



An c. 



A-l 



(27rm„T)3/2(^-i)yl3/2- 



(12) 
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By comparing equations 1)12(1 and we obtain a formal 
relation between the model parameters: 

A possible physical interpretation of this relation is the fol- 
lowing. If we assume that the freeze-out density and the 
coalescence parameter are determined by a short range in- 
teraction between nucleons and properties of formed frag- 
ments, then there is an effective 'temperature' Ta which 
characterizes the produced coalescent fragments. In the 
case of saturation of the binding energy, i.e. Ba ~ A, the 
effective temperatures do not differ much for fragments 
with different A. That resembles the chemical tempera- 
ture, when the relative probabilities of different fragments 
are determined by this temperature. One can see that the 
isotope temperature defined by equation corresponds 
exactly to this temperature. Moreover, one can also use 
relation H13|) in another way, namely, one can determine 
the coalescence parameters po for different fragments from 
the experimentally obtained chemical temperatures .'Ui . 
Therefore, if statistical models are applied to interpret 
processes with strong dynamical features, we should take 
into account that the fragment velocities may be deter- 
mined by the kinetic temperature (or by the initial dy- 
namics) but not by the chemical temperature. 

3.4 Radial flow 

After a central collision the collective expansion is as- 
sumed to change only the velocity of the fragments taken 
into account by a flow velocity profile v/(r) = (r/i?) • vq 
proportional to the position r from the center of the disas- 
sembling equilibrated source ^J. Within this scenario the 
flow velocity is superimposed onto the stochastic motion of 
the generated fragments and only the stochastic thermal 
part is responsible for fragment production. The radial 
flow is supposed to change the fragment velocities but not 
the fragment yields. This ansatz comes from the hydrody- 
namical picture, however, in the case of nuclear multifrag- 
mentation there are theoretical arguments for using this 
approach at Efiow < 3 A-MeV [T1IS7| . In this paper we 
extend it also to higher flow energies in order to find lim- 
its of this approach by analyzing the experimental data. 
The adequacy of the above ansatz is formally supported 
by statistical-like properties of fragments produced in a 
dynamical process such as coalescence. The lattice model 
calculations 38 show also that the fiow may not influence 
statistical fragment formation. 

This hypothesis of decoupling thermal and collective 
motions is sufficient for a reasonable reproduction of ex- 
perimental data PIQIEIEI ■ By introducing the phenomeno- 
logical radial flow profile we fit the IMF velocities, how- 
ever, we do not explain the velocities. Nevertheless, since 
we conserve the total energy and momentum in the sys- 
tem, this receipt allows to simulate individual multifrag- 
mentation events and to compare them directly with ex- 
perimental data. 



The described modified statistical model does not pre- 
tend to be a complete substitute for dynamical calcula- 
tions but it should be considered as a effective tool for a 
primary analysis of experimental data. Then, the observ- 
ables being properly reproduced, the physical meaning of 
the fitted parameters can be interpreted. The knowledge 
of these parameters (in particular, the temperatures) is 
supposed to be important for applications of thermal de- 
scriptions for many nuclear processes. The code SMMFC 
allows us to perform calculations with high statistics for 
large systems, as Au on Au, with reasonable computing 
time. The code produces event distributions directly re- 
lated to the observables, e.g. the multiplicity of an event 
as well as the charge, the mass, the kinetic energy and the 
polar and azimuthal angles. In this way, the generation 
of single events allow us to process the calculated quan- 
tities in a way like the experimental data samples. Our 
implemented filter permits to study the influence of the 
detector geometry and resolution, the Z-dependent regis- 
tration thresholds and chosen cut conditions. 



4 Adjustment of the model parameters 

The following four parameters of the model (see section 3) 

E* , Arei , Efiow and P-^ 

determine the fragment production. The influence of the 
break-up density g is discussed below. A scanning over 
the complete parameter space, used in [3|S], is too time- 
consuming in our case so that we disentangle this coupled 
parameter set by finding specific sensitivities of the pa- 
rameters to certain observables. We refer in the following 
only to the key observables which are necessary to fix the 
model parameters. 



4.1 Thermal excitation energy 

Multiplicity distributions of charged particles from central 
Au-|-Au collisions (ref. 13 ) are plotted in Fig. 3. In addi- 
tion, our results obtained from a subset of data selected by 
the centrahty criterion ERAT > 0.7 and Oc.m. > 25° are 
also shown in this figure. Since the experimental distribu- 
tions have a nearly exponential shape dN/dZcx exp{—aZ) 
over the whole range of measured charges, they can be fit 
by an exponential form. We included into the flts only the 
part of the charge distributions with Z > 3 to exclude the 
influence of LCP on the parameter a. In ref. |2| also the 
Z=4 points were excluded from the fits performed within 
3 < Z < 10. Our fit parameters a(3 < Z < 6) obtained 
from the mentioned subset are somewhat larger than the 
corresponding one's of ref. |13j . Since the errors Aa from 
each fit are small compared to possible disturbance of the 
charge distribution by evaporative processes and sequen- 
tial decays, we used the difference between the fit values 
given in ref. J3] and a(3 < Z < 6) as a measure of possi- 
ble deviations. 
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Fig. 3. Measured and simulated multiplicities N of nuclear 
charges. Black dots: data from ref. 113) . Open triangles: subset 
of data analysed by the selection criterion ERAT > 0.7 and 
Oc.M. > 25° within the plastic wall acceptance, Open circles: 
filtered model calculations with input parameters from table 
1. The calculated Z=l,2 multiplicities contain contents from 
both the preequilibrium and thermalized sources. 



The parameters a derived from the data are nearly 
proportional to the available center-of-mass energy. 0th- 
erwords, the calculated charge distributions show that the 
steepness parameter a{3 < Z < 6) increases also linearly 
with the excitation energy of the thermal source. In order 
to find a quantitative relation, numerous charge distribu- 
tions were simulated at various excitation energies. From 
a linear fit within the limits 3<Z<6 we found the relation 
a = —0.571 + 0.1508 • E* within the considered ranges of 
the source parameters. This result is nearly independent 
on the source size Arei , freeze-out density g and the radial 
flow Fifiow Steepness parameters a obtained from calcu- 
lated charge yields for the full phase space deviate only 
by a few percents from those which were processed by the 
filter routine implemented in SMMFC. 

The procedure of determining the parameters at 250 
A-MeV is illustrated in Fig. 4. The values a=0.91 and 
a=1.15 IE! found from fits to the data determine the 
interval of the corresponding excitation energy E* (upper 
panel of Fig. 4). 



4.2 Size of the thermal source 

The interval of excitation energies being fixed, we are go- 
ing to estimate the relative source size Arei by means of 
the IMF multiplicity which depends in terms of SMMFC 
on both Arei and E*. The lower panel of Fig. 4. demon- 
strates how Arei is estimated from the overlap of the ex- 
perimental IMF multiplicity jJS] with the calculated one's. 
Arei = 0.50 ± 0.05 at 250 A-MeV is consistent with both 
the multiplicity data and the estimated limits of E* . The 
expected sizes of the equilibrated source at 150 A-MeV 
and 400 A-MeV have been obtained analogously and the 
results are given in table 1. These parameters are consis- 
tent with those extracted in ref. jS] at somewhat lower 
incident energies. The obtained trend of the decreasing 
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Fig. 4. Determination of the input parameters of the model. 
Upper part: steepness parameter a vs. excitation energy E*. 
Dashed horizontal lines: data taken from I13| and this work. 
SMMFC calculations are denoted by symbols. Full squares: 
E/iou) =0.55 ■ Ec.A/., A,.e;=0.3. Black triangles: E/;ojd =0.55 
• Ec.Af., Arei =0.5, Open squares: E/iou, =0.37 ■ Ec.m., 
Ar.ei=0.3. Lower part: IMF multiplicity (3< Z < Q) vs. ex- 
citation energy. Hatched area: data from ref.|13|. The symbols 
are SMMFC calculations for three sets of Arei at fixed param- 
eters E* =11.0 MeV and E/;o„ =0.55 • Ec.a/.. The estimated 
range of E* is between the dotted vertical lines. 



size of the thermal source with the beam energy is sup- 
ported by dynamical calculations. For example, the anal- 
ysis of the FOPI data with the hybrid model BUU+SMM 
of ref. 9 has shown that in central Au-t-Au collisions the 
fraction of thermalized matter drops from Arei=0-^S at 
150 A-MeV to Arei=0.3 at 250 A-MeV, respectively. 
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Table 1. Properties of the thermal sources. Deviations of the 
parameters given in parentheses are admissible to reproduce 
the experimental data within the error bars.E/jou, is used in 
the model calculations. The parameters for 100 A-MeV and 1 
A-GeV are extrapolated one's. 





E* 


Tth 




Eflow 


(A-MeV ) 


(A-MeV) 


(MeV) 




(A-MeV) 


100 


~9.0 


~7.3 


~0.81 


10.5(1.5) 


150 


10.0(1.1) 


7.8(0.4) 


0.72(0.11) 


20.(3.0) 


250 


10.8(0.8) 


8.4(1.1) 


0.50(0.05) 


32.(6.) 


400 


15.0(1.4) 


11.5(2.4) 


0.50(0.13) 


56.8(6.) 


1050 


24.0(2.0) 


21.2(5.0) 


0.26(0.10) 


81.(10.) 



Table 2. Reproduction of IMF multiphcities. 



E beam 


{IMF) 


{IMF) 


(A-MeV) 


SMMFC(47r) 


data(47r) 


150 


11.61 


10.35±0.06 13 


, 12.6 I14J 


250 


8.42 


7.58±0.04 13 


,8.3 14 


400 


5.50 


5.16±0.03 il3| 


■ 5.7 jl4j 



4.3 Radial flow 

As start parameter we used results from previous analyses 
|i;-iLl l5"17 39 which were slightly varied to get an optimum 
reproduction of the fragment's mean kinetic energies. The 
values Eflow and the allowed spread given in table 1 pro- 
vide a satisfactory agreement between data and calcula- 
tion. These values are in agreement with other analyses of 
the data gOj. 

4.4 Energy release by plons 

The quantity P^^ was estimated (i) from pion multiplicities 
predicted by Vlassov-Uehling-Uhlenbeck transport model 
calculations in dependence of the beam energy at impact 
parameters of b=3 fm j^T], (ii) from the Z\-resonance pro- 
duction probability in central Au-t-Au collisions 021 : (iii) 
from the pion-to-proton ratios at 1.05 A-GeV 031 and 
(iiii) the data in ref. 05. From all those references one 
can conclude that pion production at 1.05 A-GeV carries 
away ~10 % of the available CM. energy. For the lower 
considered energies the pion contribution is neglegible. 

4.5 Break-up density 

In the model SMM/SMMFC the free volume influencing 
the translational entropy of partitions is not determined 
by the total volume (or density) of the system, though 
this assumption is adopted in some other statistical anal- 
yses The free volume reflects dynamics of fragment 
formation and it depends on the fragment multiplicity 1 . 
This ansatz is important for a good reproduction of ex- 
perimental data 10 . In the model the density influences 
directly only the Coulomb interaction in the system. The 



Table 3. Calculated multiplicities of LCP's in central Au+Au 
collisions using the scenario (ii) and the parameters from table 
1. Typical uncertainties in the resulting sum are ± 5. 





E(A-MeV) 


pre 


equ 


sum 


data '13 


Z=l 


150 


32.2 


32.4 


64.6 


61.84(0.58) 


Z=l 


250 


52.2 


30.8 


83.0 


75.82(0.62) 


Z=l 


400 


60.0 


34.8 


94.8 


92.04(0.62) 


Z=2 


150 


5.3 


15.6 


20.9 


26.76(0.36) 


Z=2 


250 


8.9 


14.5 


23.4 


27.27(0.36) 


Z=2 


400 


9.3 


14.3 


23.6 


24.16(0.30) 



Table 4. Same as table 3, but calculations without coales- 
cence. 





E(A-MeV) 


pre 


equ 


sum 


data 13; 


Z=l 


150 


40.6 


32.4 


73.0 


61.84(0.58) 


Z^l 


250 


69.0 


30.8 


99.8 


75.82(0.62) 


z=\ 


400 


77.8 


34.8 


112.6 


92.04(0.62) 


Z=2 


150 


1.4 


15.6 


17.0 


26.76(0.36) 


Z^2 


250 


0.6 


14.5 


15.1 


27.27(0.36) 


Z^2 


400 


.1 


14.3 


14.4 


24.16(0.30) 



resulting excitation energies of the equilibrated source are 
large compared to the Coulomb energies so that only mi- 
nor changes in the fragment yields are expected if the 
freeze-out density is changed. This was confirmed by cor- 
responding calculations within < q/ Qq < 1/3, where 
is the normal nuclear density. Our final calculations 
were performed with a freeze-out density of p=l/6 -qq. 



5 Comparison of experimental data with 
model calculations and discussion 

5.1 Fragment charge distributions and multiplicities 

Since our input parameters for SMMFC have been de- 
rived only from the IMF observables it is important to 
proof to what extent the model reproduces also the LCP 
multiplicities. Figure 3 and table 2 show that the model 
is able to reproduce the charge distributions from LCP's 
to IMF's in shape as well as in absolute scale provided 
that preequilibrium nucleons and the coalescenced parti- 
cles are taken into account. On the contrary, calculations 
performed without coalescence show a clear underestima- 
tion of the Z=2 multiplicity as seen by comparison of ta- 
bles 3 and 4 based on the data taken from Ref. |13| . 

5.2 Isotope Temperatures Tiso 

In order to describe the experimental findings shown in 
Fig. 1 we performed model calculations by using the pa- 
rameters evaluated in section 4 (see table 1). The calcula- 
tion at 100 A-MeV was performed with extrapolated pa- 
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Fig. 5. Temperatures as a function of the beam energy. Dots: 
isotope temperatures obtained from the yield ratios d/t and 
'^He/'^He (integrated over the particles kin. energies') |15ll6l . 
open dots: the same but in the velocity range limited by u/c « 
0.3, black triangles: the same, taken from 46 , black asterisk: 
temperature from ref. |48|, hatched area: limits of the micro- 
canonical temperature calculated with SMMFC using the input 
parameters of table 1. 



rameters given in table 1, which are consistent with pa- 
rameters extracted in ref. 3 . The generated kinetic en- 
ergy distributions of d, t, ^He and ^He were filtered by 
the cuts set in the experiment ^^l- Then the calculated 
LCP spectra were transformed into velocity distributions 
like the experimental data. The temperatures Tn-He ob- 
tained from that reproduce almost quantitatively the data 
(see Fig. 1.) 

In the calculations the increase of Tiso with velocity 
is mainly caused by a combined effect of the thermal and 
coalescence mechanisms, since they favour production of 
LCP's with different energies [221 ■ The temperature in this 
range of kinetic energies is also sensitive to the radial flow: 
the simulations undershoot the data if too less radial flow 
is assumed and overshoot them for too much flow. The 
parameters within the limits given in table 1 allow a rea- 
sonable reproduction of the data. 

The isotope temperatures T^^o obtained from energy- 
integrated yields, as well as from the yields in the limited 
range of fragment velocities (see section 2) are presented 
in Fig. 5. This figure shows also results of the ALADIN 
collaboration obtained for central Au-|-Au collisions in the 
energy range from 50 A-MeV to 200 A-MeV as well as 
data at 35 A-MeV US], which match also the found trend. 
Figure 5 shows also the microcanonical temperature Tth 
of the thermal source calculated with the code SMMFC. 
As discussed in this temperature is slightly different 
from Th-Hc^ however it shows clearly the same behaviour 
with increasing beam or excitation energy. Therefore, the 
isotope temperature at low beam energies can be used to 
deduce the temperature of the thermal source j45| . 



5.3 Kinetic energy spectra 

In this section we compare measured kinetic energy distri- 
butions of fragments produced at 250 A-MeV and 1 A-GeV 
with corresponding calculations. Calculations with only 
the thermal source fail to reproduce the tails of the spec- 
tral shapes of Z=l and Z=2 particles at 250 A-MeV (see 
Fig. 6). A consideration of the preequilibrium contribution 
improves the calculated shape. The displayed two data 
sets of Z—1 distributions demonstrate that the spectral 
shape is rather sensitive to the criterion how central events 
are selected. The calculated SMMFC distribution is close 
to the event sample obtained by the stringent combined 
criterion of ERAT and directivity (see ref. J3j). 

The disagreement of the tails in the proton spectra 
might be due to deviations of the initial distributions 
of protons from a one-particle equilibration. Neither the 
spectral shape nor the multiplicity of light clusters can be 
reproduced without consideration of coalescence. 

Generally, the presence of the radial flow affects the 
energy spectra and imitates a high temperature. As sup- 
posed, the radial velocities of the fragments depend on 
their positions in the freeze-out volume. This leads to 
an additional difference between velocities of fragments 
which is not connected with their thermal random mo- 
tion. Therefore, the deduced source temperature may be 
essentially lower than the temperature extracted from fits 
of kinetic energy distributions to the data. For example, 
the calculations shown in Fig. 6 were performed with ra- 
dial flow of Fifiow = 32 A-MeV. The spectra shapes of 
lithium clusters generated by a thermal distribution su- 
perimposed with this flow are close to the measured data. 
The Siemens- Rasmussen formula which does not take 
into account the effect of the fragment positions, repro- 
duces the shape of this distribution also, if the kinetic 
temperature is Tkin — 29 MeV. However, this tempera- 
ture is considerably larger than the thermal one. 

In Fig. 7 we compare the mean kinetic energies calcu- 
lated also for the heavier clusters with the corresponding 
data [Sni- The satisfactory agreement between data and 
calculation is consistent with the result of ref. ^21 where 
it was found that 34.0 ± 3.9 A-MeV of the available CM. 
energy is stored into the radial flow. 

Next we applied SMMFC to analyse data obtained 
in central Au-|-Au collisions around 1 A-GeV |16lll7j . At 
this incident energy mostly LCP's are produced and the 
influence of the thermal source on the fragment produc- 
tion is very limited. The model input parameters E* and 
Arei (given in table 1) were estimated from the values at 
low Ebearn by Straightforward extrapolation. The energy 
stored in radial flow was taken from ref. jl?) . A slight vari- 
ation of the parameters (±12%) was allowed to find op- 
tional agreement with the data. Proceeding on this input, 
the model calculations ^ give a quite reasonable descrip- 
tion of the proton distribution (Fig. 8), demonstrating the 

no filter was applied to the calculated events since it was 
emphasized f ref. 1161 ') that 'the TPC... allows measurements of 
spectra up to angles of 90° in the center of mass with no low-pr 
cut for central events'. 
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Fig. 6. Kinetic energy spectra of ejectiles with charges 
Z=l, 2, 3 at 250 A-MeV. FuU dots: experimental data for the 
cut 0c. M. > 25° and ERAT > 0.7, fuU triangles: data from 
ref.|13|. Dashed histogram: thermal spectrum, open circles: 
sum of thermal and preequilibrium parts (phase space gener- 
ation, coalescence included). Calculations and data are scaled 
among themselves to compare the spectral shapes. 
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Fig. 8. Differential kinetic energy distributions of protons at 
1.0 A-GeV. Dots: EOS data (17| . dotted line: calculated spec- 
trum of the thermal source, solid line: sum of the fast-stage and 
thermal sources normalized to the data, dashed line: fit of the 
calculated proton distribution using the Siemens-Rasmussen 
relation. 



Table 5. Multiplicities of LCP's in central Au+Au collisions 
at 1.0 A-GeV. Upper row: data |16I18| . lower row: SMMFC 
calculations. 



protons 


deuterons 


tritons 


^He 




78.82 
±1.48 


33.34 
±0.98 


9.56 
±0.48 


6.06 
±0.4 


2.40 
±0.26 


84.7 


35.9 


9.6 


9.7 


3.6 



I ^ ^ ^ ^ ^ ^ 1 

1 2 3 4 5 6 7 

NUCLEAR CHARGE Z 

Fig. 7. Mean kinetic energies per nucleon at 250 A-MeV. Dots: 
data from ref. |39| for the azimuthal angle (p—180°. Hatched 
area : calculations with 32< Efiow < 34 A-MeV. 



high degree of one-particle equihbration reached in this 
central event sample. However, protons from the complete 
thermalized source contribute only with a very small frac- 
tion (see the dotted line in Fig. 8). In tables 5 and 6 we 
compare experimental and calculated multiplicities and 
average kinetic energies of LCP's. A sufficient agreement 
between data and calculations could be achieved. The data 
given in table 5 were used to calculate the temperature 
Th-Hc = 18.2±5.9 MeV. The corresponding isotope tem- 
perature calculated with SMMFC amounts to T^so — 18.0 
MeV. However, these temperatures are considerably lower 
than the kinetic temperatures determined by the slope of 
the energy spectra. 
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Fig. 9. Energy distributions of hydrogen and hehum isotopes. 
Black dots: data taken from ref. |17| . Solid histograms: sum of 
the fast-stage and thermal sources. Input Efiom =90 A-MeV, 
Are;=0.26, P,r —0.10. Dashed histograms: the same input, but 
without coalescence. Calculations are normalized to the data. 



Table 6. Average kinetic energies < Ekin > in central Au+Au 
collisions at 1.15 A-GeV. Upper row: data |17|. lower row: 
SMMFC calculations. 



protons 


deuterons 


tritons 


^He 


^He 


203±3 
MeV 


263±3 
MeV 


322±12 
MeV 


328±12 
MeV 


359±20 
MeV 


233±11 
MeV 


278±18 
MeV 


331±12 
MeV 


307±15 
MeV 


317±15 
MeV 



In order to illustrate this point we compare our results 
with the fit of the data performed in |16II17| based on the 
blast scenario of Siemens and Rasmussen Q^l- In Fig- 8 
we show also the corresponding blast model fit to the pro- 
ton spectrum equivalent to the high kinetic temperature 
of Tfci„=81 MeV 22]. Nearly the same temperature repro- 
duces also the d, t,^ He and '^He distributions nT|. The 
obvious difference between kinetic and isotope tempera- 
tures suggests that each of both temperature evaluations 
is related to different processes. 

In Fig. 9. we compare the differential LCP spectra 
taken from ref. (HI with SMMFC calculat ions using iden- 
tical input parameters. A reasonable simultaneous repro- 
duction of the spectral shapes of the above species can 
be achieved, if coalescence is taken into account. Results 



without coalescence (dashed histograms) deliver shapes 
quite different from the data. The shapes of LCP energy 
distributions calculated by SMMFC including coalescence 
are similar to the corresponding one's obtained by the 
blast-model fits. From that one can conclude that the ki- 
netic temperature, evaluated from such fits, characterizes 
the initial distribution of nucleons but it is different from 
the isotope temperature characterizing the chemical com- 
position of the produced fragments. 

Using the procedure described in section 2 we obtained 
the isotope temperature Th-Hb as function of the ra- 
dial velocities. However, contrary to the increasing trend 
shown in Fig. 1, at 1 A-GeV beam energy the temperature 
Th-Hc does not increase with increasing velocity. Prob- 
ably, this is a consequence of the coalescence mechanism 
which is mainly responsible for the LCP production (see 
also discussion in section 5.4). 



5.4 Correlation of the radial flow with the chemical 
temperature 

By comparing Fig. 5 with table 1 one can see that both 
the isotope temperature, characterizing the produced frag- 
ments, and the radial flow, reflecting the dynamics of 
the process, increase with the beam energy. We suggest 
that the correlation between these parameters may pro- 
vide complementary information on mechanisms of frag- 
ment production. 

As pointed out, there are two contributions to the iso- 
tope temperature given by the SMMFC. The first one is 
related to the fragment production in the thermal source, 
which dominates at low flow energies. Experimentally, this 
temperature could be identified approximately at small 
flow velocities, as seen in Fig. 1. In the following we de- 
pict Th-Hc extracted at fragment velocities v <5 cm/ns 
as the 'thermal' isotope temperature. Contrary to this 
temperature we call isotope temperatures obtained from 
energy-integrated yields 'total' ones. Figure 10 presents a 
phenomenological relation between the isotope tempera- 
ture and the radial flow energy of the thermal sources. 
The 'thermal' isotope temperatures, shown in Fig. 10 as 
open circles, increase with increasing radial flow. With re- 
gard to other experimental data we note that the isotope 
temperature corresponding to central Au on Au collisions 
at 35 A-MeV 0Hj should be considered also as the 'ther- 
mal' one, since the dominating thermal source includes 
the small radial flow energy of < 1 A -MeV 2,. 

The second contribution to the isotope temperature 
comes from the coalescence mechanism and it provides 
higher values of the temperature. In order to compare both 
contributions, we show in Fig. 10 the total isotope tem- 
peratures obtained from the same data. In addition, we 
included into this figure also 'total' isotope temperatures 
Th-Hc obtained in central Au-|-Au collisions at 50, 100, 
150 and 200 A-MeV taken from ref. ^6}. The correspond- 
ing radial flow energies were estimated by interpolation 
of the values given in table 1 and they are close to other 
results Umi. 
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The found difTerence between the 'total' and 'ther- 
mal' isotope temperatures indicates the existence of the 
two contributions. At small flow energies the difference 
increases with increasing flow, because the contribution of 
the coalescence becomes larger. This difference is supposed 
to decrease at large flow energies, since the the contri- 
bution of the completely thermalized source disappears. 
At 1 A-GeV beam energy, IMF are scarcely produced. 
Therefore, there are large uncertainties in both the radial 
flow and the isotope temperature contributions obtained 
for the thermal source. However, within the error bars no 
difference between the 'total' and the 'thermal' chemical 
temperatures was found. That is consistent with the dis- 
appearance of the thermalized source. The comparison of 
the 'thermal' and the 'total' temperatures shown in Fig. 10 
suggests that the transition between the described mech- 
anisms accomplishes rather smoothly since the tempera- 
tures are not very different. 

Figure 10 shows also data obtained from another type 
of reactions, namely break-ups of Au nuclei after periph- 
eral collisions These data correspond to small flow 
energies and they match the general trend indicating that 
the fragment formation seems to correlate with the ap- 
pearance of the radial flow independently of its origin. 

From the obtained results we suggest the following evo- 
lution of the fragment production process with increasing 
radial flow. As long as the flow is small, the fragments 
are mainly produced in a completely equilibrated source. 
However, the energy available for the thermal population 
of the phase space does not include the flow energy. In- 
creasing radial flow influences the fragment formation in 
a twofold way: (i) it increases the velocities of nucleons 
forming a fragment, and (ii) it restricts the phase space 
population by cutting many-particle correlations between 
nucleons. This leads to the production of small fragments, 
which could be effectively described as increasing temper- 
ature of the thermal source. At larger radial flow, only 
one-particle correlations remain, which correspond to the 
coalescence mechanism. Nevertheless, properties of com- 
posite particles produced by coalescence resemble features 
of statistical processes. In particular, one may introduce 
a common temperature for any of the fragments which 
characterizes their formation probabilities and which is 
an ingredient of the chemical equilibrium. Therefore, one 
can formally proceed to treat the fragment formation sta- 
tistically, but taking into account the specific relations be- 
tween the new 'statistical' parameters caused by the dy- 
namics of the process (see, e.g., relation (jlHIl '). Some rela- 
tions of the standard thermodynamics, such as the equiv- 
alence of the chemical and the kinetic temperatures, are 
not valid in this case. However, one can deduce relations 
between statistical and dynamical parameters, like that 
shown in Fig. 10. 

6 Conclusions 

We have developed an extended version of the Statistical 
Multifragmentation Model (SMMFC) aimed at the analy- 
sis of experimental data related to LCP and IMF produc- 
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Fig. 10. Isotope temperatures versus radial flow. Open cir- 
cles: thermal temperatures estimated at low fragment veloci- 
ties. Open triangle: thermal temperature taken from 148 1 . Full 
symbols correspond to temperatures obtained from energy- 
integrated isotope yields. Dots: FOPI data |15) . Triangles: AL- 
ADIN data gg. Square: EOS data |16I17| . Full asterisks: pe- 
ripheral collisions from |10|. Open asterisks: peripheral colli- 
sions from |11I47| . Dashed line: Calculated isotope tempera- 
ture of the thermal source alone. Solid line: isotope tempera- 
ture Tu-He calculated by SMMFC including coalescence. 



tion in central heavy-ion collisions. In the present studies 
the model was applied to analyse Au-f Au data obtained 
by the FOPI and EOS collaborations in the energy range 
from 100 to 400 and at 1000 A-MeV, respectively. 

The statistical analysis includes the following physical 
processes of the reactions. An essential part of nucleons of 
the whole system is released during the fast (dynamical) 
stage and the remaining matter constitutes an equilibrium 
source. The share of the complete thermalized source Arei 
decreases with increasing beam energy. This result is in 
agreement with analyses of experimental data carried out 
in refs. [2111113 ■ It was found that the thermal sources have 
temperatures which are considerably lower than expected 
from the kinetic energies of the produced fragments. Both 
the nucleon and the LCP emissions in the fast stage are 
suppossed to be the reason of the reduced equilibrated 
energy at the freeze-out. Moreover, an appreciable part 
of the available energy of the equilibrated source is con- 
verted into collective motion, e.g. radial flow. Though the 
analysis remains phenomenologically in part, the found 
regularities seem to be reliable since they are supported 
also by analyses of other data. Therefore, the results may 
be used for interpolation and qualitative estimations of 
parameters of thermal sources. 

According to our findings, a considerable growth of the 
flow energy is accompanied by a very moderate increasing 
of the thermal temperature in these central collisions. In 
our opinion, the most reliable determination of the tem- 
perature in this case should be achieved by the isotope 
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thermometer which can be directly related to the thermal 
source. As we have shown, it is possible to apply this ther- 
mometer for testing the dependence of the temperature on 
the fragment velocities, which can be used to identify dif- 
ferent mechanisms of fragment production. 

Light clusters can also occur as result of a dynamical 
process which involves secondary interaction of the fast 
nuclcons. We have shown that the coalescence mechanism 
is responsible for the production of light fragments and 
its contribution dominates at larger beam energies. Coa- 
lescence is caused by the short range attractive interac- 
tion between the nuclcons, and, imdcr some conditions, it 
is consistent with the chemical equilibrium. Nevertheless, 
the isotope temperatures obtained from yields of coales- 
cent fragments remain very moderate in comparison with 
their kinetic energies (i.e. with their kinetic temperatures). 
The analysis of experimental LCP energy distributions 
leads also to this conclusion. Therefore, one can speculate 
that the chemical equilibrium is attained in such exploding 
systems. However, it is different from the one-particle ki- 
netical equilibration. We have demonstrated that the iso- 
tope tempcratme is correlated to the radial flow in central 
collisions. This presumed relation between isotope temper- 
ature and radial flow may be used to estimate chemical 
temperatures in different explosive processes. 
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